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Abstract 

We study the single spin asymmetry (SSA) for the pion production in semi-inclusive 
^ . deep inelastic scattering, ep^ — > enX, in the framework of the collinear factorization. We 
derive the complete cross section formula associated with the twist-3 quark-gluon correlation 
functions for the transversely polarized nucleon, including all types of pole (hard-pole, soft- 
fermion-pole and soft-gluon-pole) contributions which produce the strong interaction phase 
necessary for SSA. We prove that the partonic hard part from each pole contribution satisfies 
certain constraints from Ward identities for color gauge invariance. We demonstrate that 
the use of these new constraints is crucial to reorganize the collinear expansion of the 
Feynman diagrams into manifestly gauge-invariant form so as to obtain the factorization 
formula for the cross section in terms of a complete set of the twist-3 distributions without 
any double counting. It also provides a simpler method for the actual calculation. We also 
present a simple estimate of SSA based on our cross section formula, using a model for the 
"soft-gluon-pole function" that represents the relevant twist-3 quark-gluon correlation, and 
compare the magnitude of the terms involving the derivative of the soft-gluon-pole function 
with that of the "non-derivative" terms. 
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1 Introduction 



Understanding of the large single transverse spin asymmetries (SSAs) observed in pp col- 
lisions in 70s and 80s [1, 2] has been a big challenge for QCD theorists. (See [3] for a 
review.) Recent data on SSA at higher energies in pp collisions at RHIC [4, 5, 6], and 
SSAs in semi-inclusive deep-inelastic scattering (SIDIS) [7, 8] have been providing even 
more exciting opportunities to clarify fundamental mechanisms of SSA as well as the rel- 
evant hadron structures. Time-reversal invariance in QCD implies that the nonzero SSA 
is caused by the interference of the amplitudes which have different phases. In addition, 
the nonzero SSA requires a helicity-flip mechanism induced by certain mass scale. In the 
framework of the factorization formula with conventional twist-2 parton distribution func- 
tions and fragmentation functions, the corresponding effects are provided through purely 
perturbative mechanism and is known to give only tiny SSA [9]. Therefore, the observed 
large SSAs require the extension of the factorization formula for high energy semi-inclusive 
reactions to incorporate novel nonperturbative functions beyond conventional twist-2 distri- 
bution/fragmentation functions, and those new functions should represent the single helicity 
flip mechanism reflecting the dynamics of chiral symmetry breaking. 

In the literature, two QCD mechanisms have been developed to explain the large SSAs, 
in conjunction with the corresponding novel nonperturbative functions. One is based on 
the use of so-called "T-odd" distribution and fragmentation functions which have explicit 
intrinsic transverse momentum k± of partons inside hadrons [10, 11]. This mechanism de- 
scribes the SSA as a leading twist effect in the region of small transverse momentum pt 
of the produced hadron. Factorization formula with the use of /cj_-dependent parton dis- 
tribution functions and fragmentation functions has been extended to SIDIS [12], applying 
the method used for e + e~ annihilation [13] and the Drell-Yan process [14], and the uni- 
versality property of those /cj_-dependent functions have been examined in great details 
[15, 16, 17, 18, 19]. Phenomenological applications of the "T-odd" functions have been also 
performed to interpret the existing data for SSAs [20]. 

The other mechanism describes the SSA as a twist-3 effect in the collinear factoriza- 
tion [21, 22], and is suited for describing SSA in the large px region. In this framework, 
SSA is connected to particular quark-gluon correlation functions on the lightcone. This 
formalism has been applied to SSA in pp collisions, such as direct photon production 
p^p — > 7X[22], pion production p^p — > ttX [23, 24, 25, 26], the hyperon polarization 
pp — > A^X [27, 28], and SIDIS ep^ — > e7rX[29]. Although the above two mechanisms 
describe SSA in different kinematic regions, it has been known that the soft-gluon-pole 
function appearing in the twist-3 mechanism is connected to a k± moment of a "T-odd" 
function [17, 30]. More recently, it has been also shown in refs. [31, 32] that the two mech- 
anisms give the identical SSA in the intermediate pt region for the Drell-Yan process as 
well as for the SIDIS, unifying the two mechanisms. 

In our previous paper [29], we studied the SSA in SIDIS, ep^ — > enX, applying the twist- 
3 mechanism. The purpose of that study was to derive the cross section formula which is 
valid in the region of the large transverse momentum of the pion, i.e. px ~ Q ^> Aqcd 
with Q the virtuality of the exchanged virtual photon, and to present an estimate of SSA in 
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SIDIS, in comparison with the similar twist-3 mechanism that reproduces the asymmetry 
A N observed for pp^ — > nX. There we identified two kinds of terms in the twist-3 cross 
section, contributing to SSA in SIDIS: 

(A) G F (x 1 ,x 2 )®D(z)®a A , 

(B) 5q(x)®E F (z 1 ,z 2 )®a B , (1) 

where G F (xi,x 2 ) and E F (zi,z 2 ) are, respectively, the twist-3 distribution function for the 
transversely polarized nucleon and the twist-3 fragmentation function for the pion. 5q(x) 
is the quark transversity distribution of the nucleon, and D(z) is the usual quark/gluon 
fragmentation function for the pion. In [29] , we focussed on the "derivative term" of the soft- 
gluon-pole (SGP) contributions of (A) and (B) terms, which is associated with the derivative 
of the corresponding twist-3 distribution function or fragmentation function E F and is 
expected to be dominant in the large-x^- or large- region. Within this approximation, 
we also presented a simple estimate of the (A) and (B) contributions: We fix the relevant 
nonperturbative functions involved in the (A) and (B) contributions so that the similar 
twist-3 mechanism associated with those functions reproduces the data on A N in pp^ — > nX. 
It turned out that the (B) term is negligible compared with the (A) term for SSA in SIDIS, 
even if we require that each of the corresponding (A)- and (B)-type contributions can 
independently reproduce An in pp^ — > ttX. 

A formalism of the twist-3 mechanism for SSA in the framework of collinear factoriza- 
tion was first presented by Qiu and Sterman for the direct photon production [22] , and it 
was applied to various other processes [23, 24, 25, 26, 27, 28, 29, 31, 32]. In this formal- 
ism, a coherent gluon field A 11 generated from quark-gluon correlation inside transversely 
polarized nucleon participates into the hard scattering subprocess. A systematic proce- 
dure for extracting the associated twist-3 effect is the collinear expansion of the relevant 
Feynman diagrams in terms of the parton's transverse momenta k±. Thus we have to 
perform the collinear expansion of a set of Feynman diagrams for the partonic hard scat- 
tering with the extra gluon field A^ attached. On the other hand, gauge-invariant twist-3 
quark-gluon correlation functions appearing in the factorization formula of the cross sec- 
tion are eventually expressed in terms of either the transverse components of the covariant 
derivative, D a = d a — igA a with a =_L, or those from the gluon field strength tensor, 
F a+ = d a A + — d + A a — ig[A a , A + ], see (2), (3) below. Therefore, a crucial step for the calcu- 
lation is to reorganize the terms from the collinear expansion of the diagrams in terms of the 
gauge-invariant combination corresponding to D a or F a+ with a =_L in the twist-3 accuracy. 
In [22], the calculation was performed in the Feynman gauge, and the combination d ± A + 
resulting from the collinear expansion is identified as a part of F ±+ = d ± A + — d + A ± + . . ., 
without any assurance about the gauge invariance of the procedure nor the uniqueness 
of the obtained results. To justify this procedure, it is necessary to show that "— d + A ±v 
also appears with exactly the same coefficient as that for d ± A + . Note that, in the Feyn- 
man gauge, the matrix element of the type (p S\ipd L A + %l)\p S) is equally important as 
(p S\ipd + A- L i(j\p S), and both matrix elements have to be treated as the same order in the 
collinear expansion. This is because, even though (p S\ipA + ip\p S) ^> (p S^A^iplp S) in 
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the Feynman gauge in a frame with p + ^> p~,p ± , acting d L to the coherent gluon field 
in the former matrix element brings relative suppression compared with the latter with 
d + acting on the gluon field. Accordingly, it is necessary to identify and keep the d + A ± 
contribution to perform the calculation in a consistent way taking into account all terms in 
the collinear expansion up to the retained order, which was missed in the literature. This 
eventually allows us to establish for the first time that the corresponding factorization for- 
mula indeed holds for the twist-3 mechanism of SSA, guaranteeing its gauge invariance and 
uniqueness. Through this development, we can also address the relations among possible 
different expressions for the cross section. 

The purpose of writing this paper is twofold: First of all, we present all necessary ingre- 
dients to establish a systematic collinear-expansion formalism for the twist-3 calculation for 
SIDIS, ep^ — > eirX. We will clarify the relation between different definitions of twist-3 dis- 
tributions, and identify and classify all pole contributions that can generate the interfering 
phase for the partonic subprocess. Then we will pay particular attention to the consistency 
relations for the hard part of the cross section which are required from gauge invariance. 
Those relations express sufficient condition to allow us to reorganize the next-to-leading 
terms in the collinear expansion in terms of the combination d ± A + — d + A ± in a unique 
way, so that eventually justify the identification of d ± A + — > F ±+ taken in [22]. This is 
a nontrivial issue since the "_L" components of a four-vector are treated in a completely 
different way from its "+" component in the collinear expansion. We will next show that all 
kinds of pole contributions for ep 1 ' — > enX indeed satisfy the above-mentioned consistency 
relations, respectively. Based on this development, we derive the complete cross section 
formula corresponding to the (A) term in (1) including all pole contributions that produce 
SSA. The obtained cross section formula is valid in all kinematic region for large px pion 
production, and has a phenomenological relevance for, e.g., the eRHIC experiment [33]. 

The remainder of this paper is organized as follows: In section 2, after recalling the 
basic properties of a complete set of the twist-3 quark-gluon correlation functions for the 
transversely polarized nucleon and the kinematics of SIDIS, we present our formalism of 
the twist-3 calculation of SSA in SIDIS. We keep all the subleading contributions of twist- 
3 from all diagrams for the cross section, and show that they can be expressed in terms 
of the twist-3 gauge-invariant correlation functions introduced earlier if the corresponding 
hard part satisfies certain constraints. The twist-2 cross section formula for the unpolarized 
SIDIS is also given for completeness. In section 3, we derive the complete result for the 
factorization formula of the twist-3 cross section corresponding to the (A) term in (1). To 
this end, we shall prove that all kinds of the pole contributions indeed satisfy the constraints 
discussed in section 2 by using the Ward identities. We also provide an economic way 
of actual calculation, and discuss the relation between different expressions for the cross 
section. In section 4, we will present a numerical estimate of SSA in SIDIS, using our cross 
section formula combined with a simple model for SGP functions. We will include all terms 
from the SGP contributions, the so-called "derivative" and "non- derivative" terms of the 
SGP functions, to see the importance of the latter contribution which was not included in 
our previous study [29]. Brief summary and conclusion are presented in section 5. 
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2 Formalism for the Twsit-3 mechanism in QCD 



2.1 A complete set of Twist-3 distribution functions 

We first recall from [29] the definition and the basic properties of the twist-3 distribution 
functions for the transversely polarized nucleon. They are characterized by the explicit 
participation of the gluon field in the lightcone correlation functions for the nucleon. For 
the transversely polarized nucleon, it is known that there are two independent twist-3 
quark-gluon correlation functions [34]. Actually, two types of different definitions for those 
functions have been used in the literature, depending on whether the explicit gluon is 
introduced through the field strength tensor F alS of the gluon or the transverse components 
of the covariant derivative, D\ = <9" — %gA\. For the transversely polarized nucleon with 
momentum p^ and spin S±, they are given [22, 34] as the "F-type" functions, 

J^J ^-e^e^-^ip SJ^(0)[0,H^(HM/^ Xn]A(Xn)\p S ± ) 

= ^ 00« e apnS -G F ( Xl , x 2 ) + i^f- ( 78 y SlG F (x u x 2 ) + ..., (2) 
and the "D-type" functions, 

/ ^ A *V^ 2 -- ){p SJ^(0)[0,H^M[/^, An]^(An)|p S ± ) 

= ^ Wy e apnS -G D ( Xl ,x 2 ) + ( 75 ^) .. S1G D ( X1 , x 2 ) + ■ ■ ■ , (3) 

where ip is the quark field, the nucleon's momentum p can be regarded as light-like (p 2 = 0) 
in the twist-3 accuracy, n M is the light-like vector (n 2 = 0) with p ■ n = 1, e apnS± = 
e aXflu p x n^Sj_ u with e i23 = 1, and the spin vector satisfies Sj_ = —1, S± ■ p = S± ■ n = 0. 
[fin, An] = P exp [ig f£ dtn • A(tn)} represents the gauge-link which guarantees gauge invari- 
ance of the nonlocal lightcone operator, and "• • •" denotes twist-4 or higher-twist distribu- 
tions. These four correlation functions, Gf(xi,x 2 ), Gp(xi,x 2 ), Go(xi,x 2 ), Gd(ii,x 2 ), are 
defined as dimensionless by introducing the nucleon mass which is a natural scale for 
chiral-symmetry breaking. Each correlation function depends on the two variables X\ and 
x 2 , where X\ and x 2 — x x are the fractions of the lightcone momentum carried by the quark 
and the gluon, respectively, which are outgoing from the nucleon. From hermiticity, P- and 
T-invariance of QCD, these functions are real, and satisfy the symmetry properties 

G F (xi,x 2 ) = G F {x 2 ,xi), G F (xi,x 2 ) = -G F (x 2 ,x 1 ), (4) 

G D (x 1 ,x 2 ) = -G D (x 2 ,x 1 ), G D (x 1 ,x 2 ) = G D (x 2 ,x 1 ). (5) 

The basis of twist-3 quark-gluon correlation functions for the transversely polarized nucleon, 
defined as (2) and (3), is overcomplete. In [29], the relation between the F-type functions of 
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(2) and the F-type functions of (3) was also established using the QCD equations motion. 
They are connected by the relations 

G D (x 1 ,x 2 ) = P — - — G F (xi,x 2 ), (6) 
X\ — x 2 

G D (x 1 ,x 2 ) = 6(x 1 - x 2 )g(x 1 ) + P — - — G F (x u x 2 ), (7) 

X\ — x 2 

where P denotes the principal value, and g{x\) can be expressed in terms of the F-type 
distributions G F (x,y) and G F (x,y), and the twist-2 quark helicity distribution Aq(x) of 
the nucleon (see [29] for the detail). In the RHS of (6), the term proportional to 5(xi — x 2 ) 
vanishes due to the anti-symmetry of Gd(xi, x 2 ) as given in (5). 

It is convenient to choose the F-type functions {G F (xi, x 2 ), G F (xi,x 2 )} as a complete 
set of twist-3 quark-gluon correlation functions for the transversely polarized nucleon: From 
(6) and (7), we note that the F-type functions are more singular than the F-type functions 
at the soft gluon point x\ — x 2 , while they are proportional to each other for x\ ^ x 2 . As will 
be discussed in our analysis for ep^ — > enX below, the cross section receives contributions 
from the "soft gluon point", X\ — x 2 , of the twist-3 distribution; we assume that there is 
no singularity in the F-type functions, in particular, that G F (xi,x 2 ) is finite at x\ = x 2 
(G F (xi,Xi) = 0, due to (4)), and use the F-type functions to express those contributions 
in the cross section. Other terms which receive contributions with x\ ^ x 2 are expressed in 
terms of either F-type or F-type functions without any subtlety. 



2.2 Kinematics for ep^ —> enX 

Here we summarize the kinematics for the SIDIS, e(£)+p(p, S±) — > e(£')+ir(Ph)+X. 1 (See 
refs. [35, 36] for the detail.) We have five independent Lorentz invariants, S ep ,Xbj,Q 2 , Zf, 
and q^. The center of mass energy squared, S ep , for the initial electron and the proton is 

S ep = (p + £) 2 ~2p-e, (8) 

ignoring masses. In the twist-3 accuracy, one can set the masses to be zero for all the 
particles in the initial and the final states. The conventional DIS variables are defined in 
terms of the virtual photon momentum q = I — I 1 as 

x bj = ^' Q 2 = -q 2 = -(t-n 2 - (9) 
zp ■ q 

For the final-state pion, we introduce the scaling variable 

z f = . (10) 

p-q 



1 In this paper, we use different notation for the kinematic variables compared with our previous pa- 
per [29]. 
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Finally, we define the "transverse" component of q, which is orthogonal to both p and P h : 

p ■ p h p- p h 

q t is a space-like vector, and we denote its magnitude by 



qr = y/-q? • (12) 

To derive the cross section, we work in a frame where the 3-momenta q and p of the virtual 
photon and the initial proton are collinear, and we denote the azimuthal angle between the 
lepton plane and the hadron plane as <fi. 

For the actual calculation, it is convenient to specify the frame further. We work in the 
so-called hadron frame [35] , which is the Breit frame of the virtual photon and the initial 
proton: 

g" = (0,0,0,-Q), (13) 

Hi'°- '&)- (14) 

In this frame the outgoing pion is taken to be in the xz plane: 

Ph - 2 [ 1 + Q*> Q 'V V • (15) 

As one can see, the transverse momentum of the pion in this frame is given by PhT = z^qx- 
This is true for any frame in which the 3-momenta q and p are collinear. The lepton 
momentum in this frame can be parameterized as 

F = — (cosh ■?/>, sinh-?/> cos 0, sinh^ sin 0, — 1) , (16) 



where 

cosh^ = ^|^-l. (17) 

Q 2 

We parameterize the transverse spin vector of the initial proton S± as 

S» = (0, cos $ 5 , sin $ 5 ,0), (18) 

— * 

where $5 represents the azimuthal angle of S± measured from the hadron plane. With the 
above definition, the cross section for ep" — > eixX can be expressed in terms of S ep , Xbj, Q 2 , 
z fi Qti $ an< l ®s m the hadron frame. Note that <fi and $5 are invariant under boosts in 
the ^-direction, so that the cross section presented below is the same in any frame where q 
and p are collinear. 
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2.3 Collinear expansion and gauge invariance 

The differential cross section for ep" — > enX can be written as 

1 d 3 P d^F' p A 

<>" = 2s7„Wm WT^ 1 ' e ' )w ^"' «• Fh) ' (19) 

where 

L^(£, £') = 2(£/ u + £ V Q - g, u Q 2 (20) 

is the leptonic tensor for the unpolarized electron and W^ v is the hadronic tensor. In the 
present study we are interested in the contribution from the twist-3 distribution for the 
initial proton combined with the twist-2 unpolarized fragmentation function for the final 
pion ((A) term in (1)). Accordingly, it is straightforward to factorize the fragmentation 
function from the hadronic tensor as 

/d 7 P 
^(*)<(p,?,-r)' ( 21 ) 

3=1,9 Z Z 

where Dj(z) (j = q, g) is the quark and gluon fragmentation functions for the pion, with 
z being the momentum fraction. Since the calculational procedure is the same for the two 
terms in (21), we consider the case for the quark fragmentation in detail and omit the index 
j from w^ v below. To extract the twist-3 effect in w^ u contributing to single-spin-dependent 
cross section, one needs to analyze the diagrams shown in Fig. 1 for the hadronic tensor. 
To this end, we work in the Feynman gauge. In Fig. 1, the partons from the nucleon, which 
is represented by the lower blob, undergoes the hard interaction with the virtual photon in 
the middle blob, followed by the fragmentation into the pion as represented by the upper 
blob. Note that the fragmentation function corresponding to the upper blob has already 
been factorized as (21). Corresponding to (a) and (b) of Fig. 1, where the momenta for 
the parton lines connecting the blobs are assigned, respectively, w^ v can be written as the 
sum of the two terms: 

«V(P,?,-) = / ^-Jr\M^(k)S^(k, q A) 
z J ( Ztt } l z 



(27 



(22) 



where the superscripts "(0)" and "(1)" indicate the number of coherent gluon line connecting 
the lower and middle blobs. M^(k) and M^ a (ki,k2) represent the lower part of (a) and 
(b) of Fig. 1, respectively, and they are given as the nucleon matrix elements: 

M§\k) = J d^e^(p S ± \^(0)UO\P S ± ), (23) 
M^(k u k 2 ) = J d^J rfV M e* (fe - fcl) > S ± \^(0)9^(V)MO\P S ± ), (24) 
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(a) (b) 

Figure 1: Generic diagrams for the hadronic tensor of ep^ — ► enX, decomposed into the 
three blobs as nucleon matrix element (lower), pion matrix element (upper), and partonic 
hard scattering by the virtual photon (middle). The first two terms, (a) and (b), in the 
expansion by the number of partons connecting the middle and lower blobs are relevant to 
the twist-3 effect induced by the nucleon. 



where the spinor indices i and j are shown explicitly; S^(k,q, — ) and S^(ki, k 2 , q, - k ) 
represent the corresponding hard parts with the Lorentz indices [i and v suppressed here and 
below for simplicity, and are matrices in spinor space. Tr[- ■ •] indicates the trace over Dirac- 
spinor indices while the color trace is implicit. Since the twist-2 unpolarized fragmentation 
functions Dj(z) in (21) are chiral-even quantities, only chiral-even Dirac matrix structures 
in M^\k) and M^ a (k 1 , k 2 ) give nonvanishing contribution to Tr[- • •]. 

In the leading order in QCD perturbation theory, the hard blob for S^(k, q, — ) in Fig. 
1 (a) stands for a set of cut Feynman diagrams for the partonic Born subprocess with a hard 
gluon or quark in the final state, while the additional, coherent gluon participates in those 
partonic processes for Sjp(ki, k 2l q, in Fig. 1 (b). A standard and systematic procedure 
to identify the twist-3 contributions to w^ v of (22) utilizes the collinear expansion of the 
hard parts S^(k, q, and Sjp(hi, k 2 , q, with respect to the partonic momenta k and 
ki t 2 around the momentum p of the parent nucleon. This reads for S^(k, q, — ) 



S^(k,q,P h /z) = S^\xp,q,P h /z) + 



dS^{k,q,P h /z) 



dk c 



u^k 13 + ■ 



(25) 



k=xp 



where x = k ■ n, u af3 = g a/3 — p a n^, and "• ■ •" denote the terms of O ((k±) 2 ) which are 
twist-4 or higher. Substituting this expansion, the first term in the RHS of (22) becomes 



dx Tr 



M {0 \x)S {0 \xp,q,P h /z) 



dx Tr 



dk a 



k=xp 



,(26) 



where the integration over k and k has been performed, with oo^k 13 of (25) transformed 
into the derivative iu> a ^d^ on the quark field of (23) by partial integration, and the relevant 
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nucleoli matrix elements have now become the lightcone correlation functions: 

dX 



M§\x) = J -e iXx {p S ± \^(0)^(Xn)\p S ± ), 



ly± d I] 



2^ 



e ax ( P S ± \^(0)d a ^(\n)\pS A 



(27) 



(28) 



Here the matrix elements M^°\x) and Mo^ a (x) for the transversely polarized nucleon can 
be decomposed into several independent Dirac matrix structures based on Lorentz, P- and 
T-invariance. Among them, we only need to consider the chiral-even Dirac matrix structures 
as noted above. We find M^\x) = , y^±gT(x)/2 + . . . with the chiral-even twist-3 quark 
distribution gr(x) [37], up to irrelevant (chiral-odd or twist-4) terms, but the first term in 
(26) with the substitution M^(x) —> 'j 5 $±g T (x)/2 does not contribute to the cross section: 
The spinor trace Tr[- • •] involving 75 produces the factor % to the first term in (26) and 
the result is contracted with the real tensor of (20) to derive the contribution to (19); 
but this yields a real quantity for the cross section only when another factor % is provided 
by the hard part S^(xp,q, Phj z), which is impossible for the Born subprocess. Similarly, 
it is straightforward to see that the second term of (26) does not contribute to the cross 
section to twist-3 accuracy, noting that the second term of (26) contains the factor % and 
considering a similar decomposition of uj^M^ 13 (x) into Dirac matrix structures as in (3). 

Therefore the twist-3 cross section for ep^ — > enX arises solely from the second term of 
(22). To analyze it, we first note that the matrix element M^ CT=± of (24) is suppressed by a 
power of l/p + compared to M^ (7=+ , and gives sub leading contribution to w^ u . Accordingly, 
in the twist-3 accuracy, collinear expansion is necessary only for S+* = S^p cr /p + which 
couples with in (22). We define 



(h, k 2 , q, P h /z) = (ki, k 2 , q, P h /z)p°. 
Collinear expansion for S^(ki, k 2 , q, Ph/z) with respect to k± >2 around p gives 



(29) 



S^\ki, k 2 , q, Ph/z) = S W (x lP ,x 2 p,q,P h /z) + 



dS^(ki, k 2 , q, Ph/z) 



+ 



dSm(h,k 2 ,q,P h /z) 



dk? 



dkf 



u a pl4 + 



ki=Xip 



(30) 



hi — %iP 



where Xj = k { ■ n with % — 1, 2. Using (30) in the second term of (22), one obtains 



w 



J dx\ J dx 2 Tr 
J dxi J dx 2 Tr 



+ 



+ 



cu a M^(x 1 ,x 2 )S^{x 1 p, x 2 p, q, P h /z) 



J d Xl j 



dxoTr 



x l ,x 2 ) 



dkf 



-XiP 
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+ 



J dx\ j dx 2 Ti 



dk% 



ki=XiP 



(31) 



up to twist-3 contributions, and, similarly to (27) and (28) above, the relevant nucleon 
matrix elements have now become the correlation functions on the lightcone: 

M^(x,,x 2 ) = f^J^ e^e^-^(p S ± \i>(0)9A a (tm)i>(\n)\p S ± ), (32) 
M^(x u x 2 ) = j^-J^- e ^ e i^-^)(p sS^)9^n)n a ^{\n)\p S ± ), (33) 

M$P(x 1 ,x 2 ) = j^J^ e^e^-^ip S ± $(0) {g (d 13 A° \im)n„) ij>(\n) 

+gA a (fm)n a d /3 ilj(Xn)}\p S ± ). (34) 

The first term of (31) does not contribute to the single-spin-dependent cross section by 
the same reason as the first term of (26): Actually the former is absorbed into the latter, 
providing the 0(g) contribution of the gauge-link operator to be inserted in-between the 
quark fields in M^(x) of (27), as can be seen straightforwardly using the Ward identity 
for the coupling of the scalar polarized, coherent gluon field n ■ A into the corresponding 
hard part (see also [38]), 



S {1) (x 1 p, x 2 p,q,P h /z) 



SM(x 2 p,q,P h /z) S^(x lP ,q,P h /z) 



xi — x 2 + ie 



x 2 — Xi — ie 



(35) 



We rewrite the remaining terms of (31) as 
«V = f dxx f dx 2 TTL a f3 M^(x 1 ,x 2 )S^(x 1 p,x 2 p,q,P h /z) 



+ 



j d Xl j dx 2 Tr L° fi (M^( Xl ,x 2 ) - Mf(x b x 2 )) dS^k^P^ 
+ [ dx x [ dx^rlito^M^^x^x^ 



ki=Xip 



X 



/ dSW(h,k 2 ,q,P h /z) 



dkf 



+ 



dSV(h,k 2 ,q,P h /z) 



ki=xip 



dk% 



(36) 



and further rewrite Mq 2 ' — M^{ p in the second term on the RHS as 



(l)/3 M {1)0 
92 ~~ 1V1 dl 

M^(x 1 ,x 2 )-M^f(x 1 ,x 2 ) = Mp^(x 1 , x 2 ) + M^(xi, x 2 ), 



(37) 
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with 



xg{(d (3 A a (nnj) - (d? A p (pn))} n a ^{\n)\p S ± ), (38) 

M (1)/3 (xi,x 2 ) = /^/^ e^e"*** 2 "* 1 ^ S ± $(0)gn ff {dfA p {nn)) ^(Xn)\p S ± ) 

= -i(x 2 -x 1 )M^(x 1 ,x 2 ), (39) 

where the second equality of (39) is obtained by partial integration. One should note the 
two terms in (38) contribute in the same order with respect to the power of the large scale 
p + : Even though, schematically, in the power counting for p + 



by one power of p + , one has 



(^(d ± A + )^j) ~ (^(d+A^Tp), 



(40) 



(41) 



since d 1 - acting on the coherent gluon field causes relative power-suppression ~ l/p + com- 
pared with d + . Using the relation (37) and (39) in the second term of (36), one obtains 



w 



+ 



J dxl J 
J dxl J 



dxoTr 



dxoTr 



oj a /3 M^(x 1 ,x 2 )^(x 2 -x 1 ) 
+S^ ) (x 1 p, x 2 p,q,P h /z)} 



ki=Xip 

dsW(Kk 2 ,q,Ph/z) 



dk% 



hi — x ip 



+ I cLd I dx 2 Ti iu a pM^{ Xl ,x 2 ) 



x 



( dSW(h,k 2 ,q,P h /z) 



{ 



dkf 



+ 



dSW(h,k 2 ,q,P h /; 



(42) 



In our analysis with a single coherent gluon connecting the hard part and the soft (nucleon) 
part as in Fig. 1 (b), Mp^(x\, x 2 ) in the first term in the RHS of (42) can be identified as 
the gauge-invariant correlation function defined as (2) with the gluon field strength tensor 
(see (38)). On the other hand, the other terms in (42) involve the matrix elements that are 
not directly related to (2). Accordingly, if those terms vanish, the first term of (42) gives 
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the cross section in terms of the F-type functions. In order that this is the case, we may 
simply require that the hard parts for each of the second and third terms of (42) vanish 
separately as 

dS^(ki, k 2 , q, Ph/z) 



(x 2 - Xi, 



dk? 



+ (x ± p, x 2 p, q, P h j z) = 0, 



dS^(h,k 2 ,q,P h /, 



dkf 



+ 



dS^(h,k 2 ,q,P h /z) 



dk% 



= 0, 



(43) 



(44) 



hi — %iP 



for a =_L, up to the twist-4 corrections. These constitute the sufficient condition to obtain 
the cross section formula in terms of the gauge-invariant correlation functions. In Sec- 
tion 3, we will see that these relations are indeed satisfied by the hard part for the partonic 
subprocess relevant to SSA. 

A comment on the previous works is in order regarding the relation (43). A formalism 
using the Feynman gauge for the twist-3 calculation of SSA was first presented in [22]. In 
[22] and all the subsequent works [23, 24, 25, 26, 27, 28, 29, 31, 32], the identification of 
d ± A + — > F ±+ in the matrix element was done in the course of the collinear expansion 
ignoring other terms in F ±+ . In [39], the relation (40) was correctly recognized (see the 
argument following eq. (86) of [39]) but the relation (41) was not noticed and the d + A ± 
term in F ±+ was discarded, (see eqs. (27a) and (27b) of [39]). Clearly this procedure is not 
justified as was emphasized around (40) and (41). It is this fact that forces us to reorganize 
the first and second terms of (36) into those of (42) and require eventually (43). Therefore, 
we emphasize the condition (43) is crucial to justify the above identification, (9 ± A + — > F ±+ . 



2.4 Pole contribution to the cross section 

As has been known from the previous studies [21, 22], the single-spin-dependent cross section 
occurs as pole contributions of the internal propagators in the partonic hard cross section. 
This can be seen as follows: Since the leptonic tensor (20) for the unpolarized lepton is 
real, the cross section receives the contributions from the real part of the hadronic tensor 
W^v of (21). When (43) and (44) are satisfied, the hadronic tensor is given by the first 
term of (42) which is proportional to the factor %. Substituting (2) into this term, we 
immediately see that the trace Tr[- • ■] of the relevant Dirac matrices, combined with$f and 
27 5 #i of the two terms of (2), produces a real factor, so that the hadronic tensor receives the 
real contributions only from the imaginary part of dSW/dk%\ ki=XiP . This requires to pick 
up imaginary part from a pole of an internal propagator in dS^ /dk^ \ki=xip- 

Analyzing the cut diagrams for the hard blob in Fig. 1(b), which represents the partonic 
Born subprocess with the additional, coherent gluon attached, it is easy to see that, when 
the momenta of "external" partons are in the collinear configuration (i.e. k\^ 2 = Xi <2 p), 
some internal lines can go on shell, with the propagators reducing to one of the following 
three factors: 

(I) — _ 1 = P — ^ m5{xi - x bj ), 
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or = P h m8(x\ — Xbj), (45) 

%2 - x bj -is x 2 — x bj 

(II) — - — = P— + in5( Xl ), or — - — = P— -mS(x 2 ), (46) 

X\ — IE X\ X2 + IE X 2 

(III) = P m5( Xl - x 2 ). (47) 

X\ — X2 + is X\ — X2 

These (I), (II), and (III) produce the pole contributions called hard pole (HP) [40], soft- 
fermion pole (SFP) [22] and soft-gluon pole (SGP) [22], respectively, when integrated over 
the parton momentum fraction. For dS^/dk^ \k i= x lP of (42), the derivative of the internal 
propagators may contribute as double pole as well as single pole, but the position of the 
double pole is the same as that of (45)- (47). In Section 3 we shall calculate all pole 
contributions to the cross section. 

2.5 Unpolarized cross section 

Here we present a brief description of the calculational procedure for SIDIS. See [35, 36] 
for the detail. Following [35], we introduce the following 4 vectors which are orthogonal to 
each other [35] : 

= ^ (9" + 2x bjP n , 



F m = ^z v x p T a , 

Z" = - q -. (48) 

These vectors become T M = (1, 0, 0, 0), X» = (0, 1, 0, 0), V = (0, 0, 1, 0), Z» = (0, 0, 0, 1) in 
the hadron frame. Note that T, X and Z are vectors, while Y is an axial vector. Since W^ u 
satisfies the current conservation, q^W^ = q v W^ u = 0, W^ v can be expanded in terms of 
9 independent tensors, for which one can employ the following: 

vr = x^x v + yr, vr = + z ll z\ vr = t»x v + xt", 

Vr = X»X V - Y"Y V , V£ v = i(T»X" - X"T U ), Vg v = i{X»Y v - Y»X U ), 

= i(T^Y u - Y^T U ), V^" = T^Y" + Y»T V , = X»Y V + VX" '. (49) 

For both unpolarized twist- 2 cross section and single-spin-dependent twist-3 cross section, 
the corresponding hadronic tensor W^ u is symmetric. In addition, for both cases, pseudo- 
tensor Vs,9 do not contribute to the expansion of W^, so that W^ u can be expanded in terms 



14 



of Vi (i = 1, • • -,4). The latter fact can be seen by explicit calculation, or by inspection: 
For the single-spin-dependent case, substituting (2) into the first term of (42), the trace 
Tr[- • •] of the relevant Dirac matrices gives the factor e h ±pn oc S± ■ (p x P^) characteristic 
of SSA, which behaves as scalar under parity transformation. 

The expansion coefficients of If" for these tensors are easily obtained by using the 
inverse tensors 14 for 14 as 

w ,u = £ vr r^vrl , (so) 



k=i 



where 



2 

= -i(T^ + xt"), vr = ^pcx" - yy). (51) 

2 2 

With these V^", we can decompose the ^-dependence of the cross sections by the contraction 
with IS" in (20). Define A k (k = 1, • • • , 4) as 

A = ^ VVf , (52) 

then one obtains 

Ai = 1 + cosh 2 ■?/>, 
A = -2, 

^4 3 = — cos0sinh2^, 

Aa = cos20sinh 2 -?/>. (53) 

With this method, the unpolarized cross section of twist- 2 for ep — > enX was derived 
some time ago [41, 35, 36]. It reads 

d 5 a a l m a s A f 1 dx f 1 dz 



dx bj dQ 2 dzfdq^,d(f) SnxljS^Q 2 £r[ 

x ]Te 2 [f a (x)D a (z)al q + f g {x)D a {z)af + f a (x)D g (z)a 9 k 



-(i-(H(H)> 

where a em = e 2 /An is the QED coupling constant, f a (x) and D a (z) are, respectively, unpo- 
larized quark distribution and fragmentation functions with quark flavor a and summation 
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J2 a is over an quark and anti-quark flavors. f g (x) and D g (z) are the gluon distribution and 
fragmentation functions, respectively. In (54), we have introduced the variables 

x = ^, z = ^, (55) 

X z 

and 

x mm = x bj (1 + TZ^fi) ' Zmm = z f( 1 + T^~^) ■ (56) 

The partonic hard cross sections in (54) are defined as (C F = (iV c 2 — l)/(2N c )) 
of = 2af = 8C F xz , 

af = AC F xz-^-(Q 2 + q 2 T ) , (57) 
Qqr 

^ = x(l-x){^(^- 2 -^ + 2) + l0-l-l) , 

[ q F \x z z z xz J x z J 



af = 2af = 8x(l - x) , 

Qq T I xz 



n 9 = 2(1 U{Q 2 + 4) - £1 , (58) 



?f = 2of = 8C F x(l - z) , 

8f = -«7^1- S ).^_{g. + ^}. (59) 
For a given S^, Q 2 and q F , the kinematic constraints for x^ and are 

< x h3 < 1 , (60) 

< 2/ < 1 - afy + x hj q\jQ 2 ' 
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Consequently, qr is limited by 



Qt < Q\ (— ~ l) (- ~ lV 

\ \ x bj J \ z f J 



(62) 



In a frame where q and p are collinear, the transverse momentum, Put-, obeys 



PhT = z f q T < z f Q 



\ \Xbj 



1 



) 




) 



(63) 



3 Calculation of the single-spin-dependent cross sec- 
tion 



In Section 2, we have seen that the single-spin-dependent cross section arises from the pole 
part of the internal propagators in the partonic hard cross sections. When the conditions 
(43) and (44) for the hard part are satisfied in the Feynman gauge calculation, the color 
gauge invariance of the cross section and the consistency of the calculational procedure are 
guaranteed. In this section, we shall show that these conditions are truly satisfied by all 
types of pole contributions in the hard part, and derive the twist-3 factorization formula 
for the single-spin-dependent cross section. Because the relevant three types of poles arise 
at the different position in the phase space as in (45)-(47), the conditions (43) and (44) 
should be satisfied separately by the corresponding hard part that arises from each pole 
contribution. We first discuss in detail each pole contribution in the hard part for the quark 
fragmentation channel, in particular, those associated with the quark-gluon correlation 
function Gf for the nucleon in Section 3.1. The results for the gluon fragmentation channel 
as well as the contributions associated with the correlation function G F will be presented 
in the subsequent sections 3.2, 3.3. 

3.1 Gp contribution: quark fragmentation channel 
A. Hard pole 

The diagrams for partonic subprocesses which produce hard-pole (HP) contributions are 
shown in Fig. 2, with fc 1)2 the momenta of "external quarks" as was defined in Fig. 1(b). The 
momentum flowing through the quark line marked by a cross is k±+q, and the corresponding 
quark propagator produces the HP at (k± + q) 2 = 0, which reduces to the pole at xi = Xbj 
after the collinear expansion, k\ — * x\p (see (45)). The corresponding HP contribution 
can be evaluated as the imaginary part of the quark propagator through the distribution 
identity, 



before the collinear expansion. In the following, we assume that the diagrams with a cross 
represent the corresponding Feynman amplitude only with its pole contribution for the 
crossed propagator. 




(64) 
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Figure 2: Feynman diagrams which give rise to the HP contributions in the quark fragmen- 
tation channel, where the hard quark fragments into final-state pion and the hard gluon 
goes into unobserved final state. The cross x denotes the quark propagator which gives 
the HP contribution. Mirror diagrams also contribute. 




Figure 3: The definition of F^, as the sum of three types of diagrams which appear in 
the LHS of the cut in Fig. 2. 

Fig. 3 shows the three types of different diagrams appearing in the LHS of the cut of 
the diagrams in Fig. 2, and we define F^(ki, k 2) q, Ph/z) as the sum of the three diagrams 
in Fig. 3, where the factors for the external lines are amputated, p is the Lorentz index for 
the coherent gluon from the nucleon, and a is the Lorentz index for the hard gluon emitted 
into the final state. F^(ki, k 2 , q, Ph/z) contains the factor 5 ((hi + q) 2 ) from the second 
term of (64). Due to the tree level Ward identity shown in Fig. 4, F^(k± : k 2 , q, Ph/z) 
satisfies the relation 

(k 2 - hf u(P h /z)F^(k 1 , k 2 , q, P h /z)t\%{k 2 + q - P h /z)5 ((k 2 + q - P h / zf) = 0, (65) 
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Figure 4: Ward identity for the HP contribution defined in Fig. 3. The dashed line repre- 
sents a scalar-polarized gluon, and the parton lines marked by a bar are on shell. 




Figure 5: The definition of fp, as the sum of two types of diagrams which appear in the 
RHS of the cut in Fig. 2. 

where u(Ph/ 

z) is the spinor for the quark fragmenting into pion, a^^{k 2 ~\~ q — Ph/ z) is the 
polarization vector for the hard gluon with the momentum k 2 + q — Phi z and the physical 
polarization A, and 5 ((k 2 + q — Ph/ z) 2 ) comes from the on-shell condition for the final 
gluon. We also define the amputated amplitude fp{k 2 , q, Ph/z) for the sum of the diagrams 
shown in Fig. 5, which represent two types of diagrams in the RHS of the cut in Fig. 2. 
Combining (k±, k 2 , q, Ph/ z) and fp(k 2 ,q, Ph/z), the sum of the six diagrams in Fig. 2 
can be written as 

Sf P (k u k 2 ,q,P h /z) 

= fp(k 2 , q, Ph/z) F^ih, k 2 , q, Ph/z) V a \k 2 + q - P h /z) 2vr 5 ((k 2 + q - P h /z) 2 ) , 

(66) 

where we have used the relation Y^ sp ms u (Ph/ z)u(Ph/ z) —Ph/z, and 

P^) = £4aW4#) (67) 

A 

is the polarization tensor for the on-shell (k 2 = 0) hard gluon with the sum of A restricted 
over the physical polarizations. Note that (66) is exactly the part of S^\k\,k 2 ,q,Ph/ z) 
appearing in (42), which comes from the HP contribution. Owing to the relation (65), 
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S p HP (ki, k 2 , q, Phi 1 z) satisfies 



(k 2 -k i ySf p (k 1 ,k 2 ,q,P h /z) = 0. 



(68) 



We now take the derivative of (68) with respect to kf (i — 1,2; a =_L) and set fcj = Xip. 
One then obtains 



(x 2 - Xi) 



(x 2 - Xi) 



dS^ p (k 1 ,k 2 ,q,P h /z)p p 



dk° 2 

dS" p (h,k 2 ,q,P h /z)pP 



+ S^ P (xip, x 2 p, q, P h j z) = 0, 



/c^ — Xip 



dkl 



-S» p (x lP: x 2 p,q,P h /z) = 0. 



(69) 



(70) 



hi — XiP 



This proves that the relation (43) holds for the HP contribution. Since x\ ^ x 2 for the HP 
contribution, the relations (69) and (70) can be written as 



dS^ p (k 1 ,k 2 ,q,P h /z)p" 



dk% 



dS" p (k u k 2 ,q,Ph/z)pP 



-XiP 



dkl 



hi — Xip 



S" P (xiP, x 2 p,q,P h /z), 



X\ — x 2 



(71) 



which proves that the condition (44) is satisfied. 

Furthermore, the relation (71) provides a useful method to calculate the corresponding 

hard cross section: Instead of calculating the derivative 



dS^ p (k 1 ,k2,q,Ph/z)p p 



dk% 



for the first 



h% — XiP 



term of (42), one simply needs to calculate S„ (xip, x 2 p, q, Phj z) j [x\— x 2 ), the RHS of (71). 
Sp P (ki, k 2 , q, Ph/z)p p represents the sum of the six diagrams in Fig. 2, and, as was noticed 
in (65), (66) above, the contribution from every diagram contains the product of two delta 
functions as 5 ((ki + q) 2 ) 5 ((k 2 + q — P h /z) 2 ). Accordingly the calculation of the derivative 

' /,J *' / ' / ' '' produces the derivative of these delta functions. However, these 



=XiP 



terms with the derivative of delta functions eventually cancel among each other, since 
Sp P (xip, x 2 p, q, Ph/z)/ (xi —x 2 ) does not contain such derivatives. This explains the absence 
of the derivatives of delta functions in the HP contribution, which has been observed in the 
literature. 2 Therefore, it is much more convenient to calculate S^ p (xip, x 2 p, q, Ph/z) / (x\ — 
x 2 ) to obtain the hard cross section for the HP contribution. 

We also note that the calculation of S^ p (x\p ) x 2 p, q, Ph/z) can be further simplified. It 
is defined as (66) in terms of the physical polarization tensor (67) which can be expressed 
as 



V af5 (k) 



-9 aP + 



k a h? + h a k? 



(72) 

it, ■ i% 

2 Thc authors of [31, 32] calculated dSf? p \h, k 2 , q, P h / ' z)p p /dk%\ k =x directly to obtain the HP con- 
tributions for Drell-Yan and SIDIS, and found the cancellation among the terms with the derivative of the 
delta functions. 
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Figure 6: Same as Fig. 2, but for the SFP contributions in the quark fragmentation channel. 



where h a is an arbitrary light-like vector (h 2 = 0) satisfying h ■ k ^ 0. However, for 
S^ p (xip, x 2 p, q, Phi z) with the transverse Lorentz index a =_L, which guarantees the phys- 
ical polarization for the coherent gluon with the momentum (x 2 — x±)p in the diagrams in 
Fig. 2, one can make the replacement V af3 (x 2 p+q — Ph/z) — > — g a/3 because the contribution 
of the second term of (72) vanishes by the Ward identities owing to the gauge invariance 
for the sum of the relevant diagrams. 



It is straightforward to check that electromagnetic gauge invariance, q^w 



fii/ 



q v w 







for our hadronic tensor (42), is satisfied by the HP contributions. Substituting (71) and 
(2) into the first term of (42), the relevant hard cross section associated with G F (x 1: x 2 ) is 
proportional to Tr S^ p (xxp, x 2 p, q, Phi ' z)]f> /(x 1 — x 2 ). Using the Ward identities of QED, 

the sum of the diagrams in Fig. 2 for Tr S^ p (xip, x 2 p, q, Ph/z)$ vanishes, when a virtual 
photon vertex is contracted by q^ or q v . 

By the method described above, it is now straightforward to calculate the HP contribu- 
tion to the cross section. 



B. Soft-fermion pole 

The diagrams which give rise to the soft-fermion-pole (SFP) contributions are shown in Fig. 
6, similarly to those for the HP in Fig. 2. The internal quark propagator which produces 
the SFP is indicated by a cross in each diagram of Fig. 6. The momentum flowing through 
this line is k\ — k 2 + Ph/z — q, so that the SFP arises at {k\ — k 2 + Ph/z — q) 2 = 0. Since each 
graph accompanies the factor S ((k 2 + q — Ph/ z) 2 ) coming from the cut propagator for the 
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X 



+ 



= 



Figure 7: Ward identity for the SFP contribution, similarly to Fig. 4. 

hard gluon, this pole reduces to that at X\ — as in (46) in the collinear limit k lj2 — > £i,2P- 
The SFP contribution can be evaluated as the imaginary part oc 5 ((ki — k 2 + Ph/ Z — <l) 2 ) 
of the corresponding propagator before the collinear expansion. As in the case of the 
Ward identity (65) for the HP contribution, the relevant Ward identity for the present 
case is shown in Fig. 7. Following the same step as for the HP contribution, it is easy 
to prove that the conditions (43) and (44) are satisfied by the SFP contribution, and also 
(71) with HP — > SFP holds. Accordingly, the SFP contribution can be calculated via 
Sf FP (xip, x 2 p, q, Ph/z)/{xi — x 2 ) and V a/3 — > —g al3 , i.e., in the same economic way as 
in the HP contribution. It is also straightforward to check that electromagnetic gauge 
invariance is satisfied by the SFP contributions using the Ward identities of QED, i.e., the 
sum of the diagrams in Fig. 6 for the calculation of Tr Sf FP (xip, x 2 p, q, Ph/z)$ / (xi — x 2 ) 
vanishes when a virtual photon vertex is contracted by q^ or q v . 

The method discussed above gives the HP and SFP contributions in terms of the F- 
type functions of (2) based on (42)- (44). It is instructive to note that the HP and SFP 
contributions can be rewritten in terms of the D-type distributions of (3) using the relation 
(6). In the new expression, S± P ' SFP (xip, x 2 p, q, Ph/z) can be regarded as the hard part 
associated with the D-type distributions. 

In this connection, we note that the authors of [22], in their study on p^p — > jX, em- 
ployed a procedure of calculating Sf FP (xip, x 2 p, q, Ph/z) directly as the hard cross section 
corresponding to the D-type distributions, and eventually obtained the results which are 
formally the same as those obtained by the above mentioned rewriting of the D-type func- 
tions in terms of the D-type functions. However, the logic for this procedure in [22] is not 
based on the relations like (71), (6), (7) in contrast to our method. In fact, the authors of 
[22] tried to find the hard cross sections corresponding to the D-type and D-type distribu- 
tions independently, and relied on their claim that the D-type functions do not receive the 
SFP contribution (see Section 3.3 of [22]). But this statement cannot be consistent with 
the relation (6) which was not noticed in [22]: If this statement were true, it would mean 
that the D-type functions do not receive the SFP contribution either, since the D-type 
and D-type functions are proportional to each other as shown in (6) and (7) for x\ ^ x 2 . 
Their argument leading to the above claim was the following: The A + component of the 
coherent gluon field represents the scalar-polarized component at the relevant accuracy in 
the collinear expansion, and thus decouples from a gauge-invariant set of diagrams in the 
hard part by means of the Ward identity, so that the hard part for the D-type function, 
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Figure 8: Same as Fig. 2, but for the SGP contributions in the quark fragmentation channel. 



{dS^(h, k 2 , q, P h /z)p p /dk 2 L - dSp(h, k 2 , q, P h /z)p p /dk^}\ k =x , identified as the coeffi- 
cient of the combination d ± A + , should vanish at X\ 7^ x 2 ■ (And the hard part for the 
F-type function can survive only for the SGP contributions with x\ = x 2 , where A + does 
not represent a scalar-polarized component.) However, the correct relations from the Ward 
identity are those as given in (69) and (70), which are different from their relation. Our 
demonstration leading to (69) and (70) indicates that it is not allowed to apply the Ward 
identities naively to the "derivative of the amplitude" which appears as a next-to-leading 
order term in the collinear expansion. We emphasize that both ingredients - - (i) the 
relations (6), (7) among the nonperturbative functions, and (ii) the correct relations (69), 
(70) from the Ward identities for the partonic subprocess - - are essential to guarantee 
the consistency of the calculation. Using only one of them would lead to contradiction or 
double counting. 



C. Soft-gluon pole 

The diagrams which give rise to the soft-gluon-pole (SGP) contributions are shown in Fig. 
8, where the quark line which produces the pole is indicated by a cross. The momentum 
flowing through this line is k\ — k 2 + Ph/z, so that the SGP arises at (ki — k 2 + Ph/z) 2 = 0, 
which reduces to the pole at x\ = x 2 as in (47) for k^ 2 — > x\$p. The SGP contribution can 
be evaluated as the imaginary part oc S ((hi — k 2 + Ph/z) 2 ) of the corresponding propagator, 
i.e., taking the quark line on-shell. 

Define Sp(k±, k 2 ,q, P^j z) as the sum of the diagrams shown in Fig. 8, where the 
factors for the external lines are amputated, the index p is for the coherent gluon, and 
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the superscript L indicates that the coherent gluon is attached to the LHS of the cut. 
Likewise we can also define S^(ki,k 2 ,q, Ph/z) corresponding to the mirror diagrams of 
Fig. 8. Sp' R (ki, k 2 , q, Ph/ z) represent the SGP contributions in the corresponding partonic 
hard scattering as directed to take the quark line on-shell by the crosses in Fig. 8. It 
also has the delta function as the cut propagator for the on-shell hard gluon. Therefore 
Sp(ki, k 2 , q, Ph/ z) is proportional to the product of the two delta functions as 



Sp(h, k 2 , q, P h /z) ~ 5 ((h -k 2 + P h /zf) 5 ({k 2 + q - P h /z) 2 ) . 



(73) 



Notice that we are considering Sp{k\,k 2 , q, Ph/ z) with parton momenta before collinear 
expansion, as in the HP and SFP cases discussed above. Since the quark-gluon vertex 
attaching the coherent gluon is sandwiched by the on-shell quark lines, the tree level Ward 
identity implies 



(k 2 -k i ysZ(k 1 ,k 2 ,q,P h /z) = 0. 



(74) 



Actually, the Ward identity of this type holds for each diagram in Fig. 8. By taking the 
derivative of (74) with respect to kf (i — 1, 2; a =_L) and setting ki = Xip, one arrives at 



(x 2 -x^ 



(x 2 - Xi) 



dS L (k l ,k 2 ,q,p h /z) P p 



dk% 



dS L Jk l ,k 2 , q ,p h /z)p p 



+ S];(x lP ,x 2P ,q,P h /z) = 0, (75) 



dkl 



hi — Xip 



S^(x lP ,x 2 p,q,P h /z)=0. (76) 



Apparently the similar relations hold for S^(ki, k 2 , q, Ph/z), too. Noting that Sp{k\, k 2 , q, Ph/z) 
+S^(ki, k 2 , q, Ph/z) is the part of S^(ki, k 2 , q, Ph/z) appearing in (42), which comes from 
the SGP contribution, the above equation (75) shows that the required consistency condi- 
tion (43) is satisfied. 

So far, the logic used for deriving (75) and (76) is the same as that used for the HP and 
SFP cases. However, there is an important difference between (75), (76) and (69), (70). 
Namely, we cannot derive the relation similar to (71) using (75) and (76), because the first 
term of (75), (76) contains singularity at x\ = x 2 due to the delta function. Owing to (73), 



dS^k 1 ,k 2 ,q,P h /z)p p /dk^ 2 



hi — %iP 



consists of the terms with 5' (xi — x 2 ) 5 ((x 2 p + q — Ph/z) 2 ), 



5 (xi — x 2 ) 5' ((x 2 p + q — P h /z) 2 ) or 5 (xi — x 2 ) 5 {{x 2 p + q — P h /z) 2 ). Because of the factor 
x 2 — Xi in the first term of (75), the calculation of S^(x lP ,x 2 p,q, Ph/z) misses the terms 
in dSp(ki, k 2 , q, Ph/z)p p /dk 2 _ , which vanish due to (x\ — x 2 )8{xi — x 2 ) = 0. 3 There- 



hi — %iP 



fore it is necessary to calculate dS^iki, k 2 , q, P h / z)p p /dk^ 
contribution to the cross section correctly. 



=x,p 



directly to get the SGP 



3 This also explains why it is difficult to obtain the SGP cross section in the lightcone gauge, A + = 0. 
In this gauge, one would identify the contribution of the F-type functions as A 1 - — > iF ±n /(xi — x 2 ) and 
calculate S± R (xip, x 2 p, q, Ph/z) to get the cross section. But a naive calculation of S± R (xip, x 2 p, q, Ph/z) 
would miss some parts of the SGP cross section as was noted above. 
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q q q q 




k 2 — fa k 2 - ki 



Figure 9: An example of a pair of diagrams from Fig. 8, with the coherent gluon attached 
on the different side of the cut. 

Because of the lacking of the analogue of (71), we have to rely on direct inspection of 
the diagrams in Fig. 8 and their mirror diagrams in order to show that the condition (44) 
is realized for the SGP contribution. We first note that, when the derivative d/dkf 2 hits 
an internal propagator which depends on k\ and k 2 only through k\ — k 2 , the condition 
(44) is automatically satisfied for that contribution. For other cases, one needs to consider 
a combination of some diagrams in Fig. 8 and their mirror diagrams. As an example, 
consider a pair of the diagrams shown in Fig. 9. The momentum flowing through the quark 
line with a cross are, respectively, k\ — k 2 + Ph/z and k 2 — ki + Ph/z for Fig. 9 (a) and (b), 
and they give SGPs which have opposite signs of the residue in the collinear limit, ki — > Xip 
(i = 1, 2). On the other hand, Fig. 9 (a) and (b), respectively, contain S ((k 2 + q — Ph/z) 2 ) 
and 5 ((hi + q — P h /z) 2 ) as the cut propagator for the on-shell hard gluon. Therefore, for 
the sum of these two diagrams in Fig. 9, the derivative hitting on these delta functions 
satisfies (44) at the SGP. If the derivative hits other internal propagator such as that with 
momentum k±+q, or k 2 + q, in Fig. 9, the corresponding contribution simply cancel between 
the diagrams, because the SGPs have opposite signs of the residue between the diagrams (a) 
and (b) in Fig. 9. It is straightforward to see that the action of the derivative d/dk" 2 can 
be classified into one of these three cases for all diagrams in Fig. 8. Therefore, the condition 
(44) is also satisfied for the SGP contribution by taking into account all the graphs including 
the mirror diagrams. This is in contrast to the HP and SFP cases, where the condition (44) 
is already satisfied by summing the diagrams which have a coherent-gluon attachment on 
one side of the cut. 

Here we mention the form of the polarization tensor from the polarization sum (67) of 
the final-state hard gluon in Sp> R (ki, k 2 , q, Ph/z). In principle, V a/3 (k 2 + q — Ph/ z) with (72) 
should be used to calculate S^(ki, k 2 , q, Ph/z), and likewise for Sjf(ki, k 2 , q, Ph/z). However, 
one is allowed to make the replacement V al3 (k 2 + q — Ph/z) — > —g alS . This is formally similar 
to the HP and SFP cases, but the logic to allow the replacement is somewhat sophisticated 
compared with those cases because we have to treat the scalar polarized coherent-gluon in 
the present case. For the proof, we first note that, in calculating the hard cross section 
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associated with the quark-gluon correlation function Gf based on (42), the corresponding 
spinor trace in the first term of (42) is taken with the insertion of the Dirac structure $f of 
the first term of (2). We consider the following replacement for this ^ in the calculation: 



2x1X2 

A difference between the original hard cross section, 

d 



(77) 



dkZ 



Tr 



S^ H (k 1 ,k 2 ,q,P h /z)p^ 



(78) 



hi — %iP 



and the new one with (77), 



d 



Tr 



S^ R (k 1 ,k 2 ,q,P h /z)p? 



2xiX 2 



(79) 



hi — %iP 



occurs when the derivative d/dk^ hits the ^ 2 in (77). However, these terms cancel between 
SpP p and S R p p for (79), because the SGP appears with opposite signs of the residue between 
them. Thus one can obtain the hard cross section using (79). We note that, in (79), \f,\ 
$2) plays the role of the factor for the incoming (outgoing) external quark line with the 
momentum k\ {k 2 ). It is now easy to see that, with the on-shell condition k\ 2 = 0, 4 the 
terms proportional to (k 2 + q — Ph/z) a or (k 2 + q — Ph/z) 13 in P a/3 (/c 2 + q — Ph/z) drops in 



Tr 



S^ R (k 1 ,k2,q,P h /z)p p ^ 1 ^2 



with the aid of the tree level Ward identity. So one 
can make the replacement V al3 (k 2 + q — Phj z) — > — g a/3 in calculating (79). Since there is no 
change between (78) and (79) even with this —g a ^, one can use (78) with the polarization 
tensor — g a/3 for the final-state hard gluon. In this way, one can safely use — g a/3 instead of 
(72) for the calculation of the SGP contribution based on (42). 

Eq. (79), combined with the on-shell condition k\ 2 = 0, is also useful to prove the 



electromagnetic gauge invariance, q^w 



q u w 



0, for the SGP contribution. Using 



the Ward identities of QED, the sum of the diagrams in Fig. 8 for the calculation of (79) 
vanishes when a virtual photon vertex is contracted by q^ or q u : The sum of vertical two 
diagrams for (79) vanishes when the left virtual photon vertex is contracted with q^, because 
of the on-shell incoming quark line with the factor k*\ and the on-shell internal quark line 
with the SGP. And the sum of the horizontal two diagrams vanishes when the right virtual 
photon vertex is contracted with g M , because of the on-shell outgoing quark line with the 
factor ^2 and the on-shell quark line emerging from the fragmentation insertion with P^. 



D. Result for the factorization formula 

Following the method described above, the single-spin-dependent cross section associated 
with the twist-3 distribution for the nucleon, Gf(x,v), and the quark fragmentation func- 

4 In general, kf = 2kfk~ — kf ± ^ (i = 1,2). Up to the twist-3 accuracy, however, we can choose as 
fcr = k^ ± /(2fc l + ), so that our external quarks are taken on-shcll kf = here. 
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tion for the pion can be obtained as 



2 



dx bj dQ 2 dzfdq^d(fi %irx 2 bj Sl p Q 2 ^ 



-7lM N 



, 1 dx r 1 dz 
sm$ s / — / — 
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X 



+ ^lG a F (x bj ,x) + a^ k G a F (0,x) 



where G a F represents Gp for quark-flavor a with electric charge e a , 5 a is defined as 5 a = 1 
for quark and 5 a = — 1 for anti-quark, and the summation J2 a is over & U quark and anti- 
quarks. The factor "—nM^/A" arises as the product of the factor % in the first term of (42), 
M N /A in (2), and in associated with evaluation of the pole contributions, crj^, 6qI, o h \ 
and a F l are the hard cross sections corresponding, respectively, to the derivative term for 
the SGP contribution, non-derivative term for the SGP contribution, HP contribution, and 
SFP contributions. They are given as 



- Vg 



4:C q q T x A Vq 
Q 2 l-f^ 



T U1 



T U2 



r V 1 



2xz 



i (Q 4 

Q 2 q\ \x 2 z 2 



+ (Q 2 



2a 



Vq 
U4 



8xz, 



'us ~ 4xz——(Q +q T ), 



(81) 
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8C q q T 
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X 
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(l-x) 2 (l-z)- 



+ 



25(2 - 35) + (1 - 2x)(l + 6z 2 - 65) 



~v q _ ~v q _ 64C q q T x 2 z 
la G4 — 



a G2 = 



~Vq 
°G3 



ICq 
Q 



-X 



Q 2 1-5 

£5(5 — Ax) 
{l-x){l-z) 



(82) 
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~Vq _ 8?T ( n . , 1 \ ^(1 + ^ 2 ) 



a H2 



a 



Vq 
H3 



o, 



g V ^ + 2A^ J (i-x)(i-i)' 



(83) 



Q 



<G F z 1 



2Nj (l-x)(l-z)' 



~Vq 



a 



F2 



.Vq 
T F3 



8C q q T 

' Q 2 



8x 2 z 
l-z 



+ x(2x- 1) + 



x 2 z 2 (2x 2 -3x + l) 
(l-x) 2 {l-z) 2 



Vq _ 64C q q T x 2 z 



Q 2 l 



' Q 



x(4x — 3) — 



xz(4x — 1) 
l-z 



(84) 



„ Vq _ 8C q q T xz(l - 2xf 



a F4 = 



Q2 (l_£)(l_i)' 

where C q = — l/2N c and o{j\ in (81) is the same as (57) for the unpolarized cross section 



(54) except for the color factor. a^ k was already given in [29]. The above result for & V G \ 



~Vq 



and oQ{ also agree with those in [32]. All the other results in (82)-(84) are new. 



3.2 Gp contribution: gluon fragmentation channel 

The twist-3 distribution G F (x,y) contributes to the single-spin-dependent cross section 
also with the gluon fragmentation function for the pion, D g (z). Again, this is classified 
into the HP, SFP and SGP contributions. The diagrams corresponding to the HP and SFP 
contributions are given by the same diagrams as shown Figs. 2 and 6, respectively, except 
that the hard gluon fragments into final-state pion and the hard quark goes into unobserved 
final state. The diagrams for the SGP contributions are shown in Fig. 10. The method of 
the calculation and the proof of the conditions (43) and (44) are the same as that described 
in the previous Section 3.1. So we shall not repeat it here. The result is obtained as follows: 

d 5 a V9 alm a s y A / -tiM n \ $ f 1 (tor 1 dz_ 

dx bj dQ 2 dz f dq^d(t) %-nx 2 hj S 2 ev Q 2 ^ x k \ 4 / S ' L mm x J Zmin z 



all [x-^G a F (x,x)\+a v G 9 k G F (x,x) 



d 



vg. 
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Figure 10: Same as Fig. 2, but for the SGP contributions in the gluon fragmentation 
channel where the hard gluon fragments into final-state pion and the hard quark goes into 
unobserved final state. 



+8ad% g k G a F (x bj , x) + 8 a a v F 9 k G a F (0, x 



Vg n a 



(85) 



where each hard cross section is given by 



~ v g __ 4C 3 g T x 



-(J. 



Vg 



Q 2 l-z" Uk) 
a v u { = 2x{l-z) 



i / Q 4 | (i-g) 2 
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Q 2 qT \ x 2 z' 



,2 1 J/n2 i ^ 2( ?T 
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(86) 
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(88) 
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(89) 



g 5 



.yg _ 8C g g T (1 - 2x) 2 x 



a F4 = 



Q 2 



1 — x 



where C g = N c /2. We remind that the hard cross section for the derivative term, (86), 
was also derived in [29] and a\j 9 k is the same as (59) for the unpolarized cross section (54) 
except for the color factor. 



3.3 The Gf contribution 

Another twist-3 distribution for the nucleon, G F (x,y), also contributes to ep^ — > enX. 
Again, in principle, the corresponding cross section arises as the HP, SFP and SGP con- 
tributions. The diagrams for those contributions are the same as in the G F case for the 
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quark fragmentation channel as well as the gluon fragmentation channel, see Figs. 2, 6, 8, 
and 10. The calculation can be carried out in parallel with the case for G F) by changing 
the relevant Dirac matrix structure to calculate the first term of (42) as $e apnS± — > i^^^S^ 
(see (2)). Owing to the anti-symmetry (4), one has G F (x,x) = 0, so that there is no 



non-derivative term for the SGP contribution. In principle, 



8Gf(xi,x) 
dx\ 



can be nonzero 



XI =x 



and could contribute to the cross section as the derivative term for the SGP contribution. 
However, after summing all the diagrams for the SGP contributions, it turns out that no 
such derivative term survives. Therefore, there appears no SGP contribution, and the result 
for the single-spin-dependent cross section associated with G F (x,y) is obtained from the 
HP and SFP contributions as 



y 2 c 



dxbjdQ 2 dzfdq^d(p 8irxljS 2 p Q : 
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,Q 2 




(90) 



where the hard cross sections for the quark fragmentation contribution are given by 
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and those for the gluon fragmentation contribution are given by 
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4 Azimuthal asymmetry 

By summing all the pole contributions with Gp and Gp derived as (80), (85) and (90) in 
Section 3, one obtains the complete single-spin-dependent cross section for ep^ — > eirX from 
the (A) term in (1). Its ^-dependence can be decomposed as 



dxbjdQ 2 dzfdqpd(p 



= sin $ s (ctq + of cos(0) + of cos(20)) . 



(95) 
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This should be compared with the similar decomposition of the twist- 2 unpolarized cross 
section of (54) [41, 35]: 



d,XbjdQ 2 dzfdqj,d(fi 



(Tq + of cos(0) + o% cos(20). 



(96) 



If one uses the lepton plane as a reference plane to define the azimuthal angle of the spin 
vector of the initial proton as (fis, and that of the hadron plane as fa, as employed in [7], 
one has the relation $5 = (fis — fa and = —(fit- 5 Accordingly, the azimuthal dependence 
of the Uq term in (95) is the same as the Sivers effect (oc sin((f) h — <j> s )) and (95) contains 
the azimuthal components which are absent in the leading order Sivers effect. 

We now present a simple estimate of SSA by using the obtained formulae. So far we 
don't have any definite information on the HP, SFP and SGP components of G a F and G a F . 
In [29], we presented an estimate of SSA including only the derivative term of the SGP 
function, -^G F {x,x), which is expected to be a good approximation at large x^ or large 
Zf. For comparison with that study, we will here include all the SGP contribution, the 
derivative and non-derivative terms of the SGP function. This calculation may be relevant, 
since we expect that there are much more soft gluons in the nucleon than the gluons with 
finite momentum. 

In order to see the relative importance of the non-derivative term, we assume the same 
model for G F (x,x) as in [29]: We take into account only the contributions from u and 
d quarks with the ansatz G F (x,x) = K a f a (x) (a = u, d) [25] and the flavor-dependent 
constant K u = —K d = 0.07, which were shown to approximately reproduce the observed 
SSA of p^p — > nX [2, 4] as SGP contributions. As in [29], we calculate the azimuthal 
asymmetries normalized by the unpolarized cross section. We define ^-integrated azimuthal 
asymmetries at $ 5 = n/2 as, 



For the purpose of illustration, we choose two sets of kinematic variables: The first one is 
S ep = 300 GeV 2 , Q 2 = 20 GeV 2 , x bj = 0.1, which is close to the COMPASS kinematics. 
Another one is S ep = 2500 GeV 2 , Q 2 = 50 GeV 2 and x bj = 0.03, which is in the region of 
planned eRHIC experiment [33]. Both sets give the same coshifj in (17). In our calculation, 
we have used GRV unpolarized parton distribution f a , g (x) [42] and KKP pion fragmentation 
function D a , g (z) [44] as in [29]. 

The results for the azimuthal asymmetries (97) for the ir° production in SIDIS are 
shown in Figs. 11(a) and (b) for the COMPASS and the eRHIC kinematics, respectively. 
The asymmetries (1)at, (cos0)jv, and (cos20)jv are shown by the solid, long-dashed, and 
dot-dashed curves, respectively, at qx = 4 GeV as a function of Zf. For comparison, 
we also showed, by the short-dashed curve, the result for (1)^ which was obtained only 
with the derivative term of the SGP contribution. One sees that (1) N is the order of 
2 % for the COMPASS kinematics and is the order of 0.5 % for the eRHIC kinematics, 

5 In [29], we made a sign mistake in this relation. 
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Figure 11: The calculated azimuthal asymmetries (97) with all the SGP contributions for 
(a) the COMPASS kinematics and (b) the eRHIC kinematics. For (1)n, we have also shown 
the result with only the derivative term of the SGP contribution. 

while (cos 0) at and (cos20)tv are negligible for both kinematics. One also sees, comparing 
the solid and short-dashed curves, that the effect of the non-derivative terms of the SGP 
contribution is important in the small and intermediate Zf region, while the derivative 
term dominates in the large Zf region. In the large Zf region, where the large x region 
of the distribution function becomes important (see (56), (80), (85)), the derivative term 
(d/dx)G F (x, x) together with our ansatz G%(x, x) = K a f a (x) causes a steep rising behavior 
of the asymmetry. If one had used a model for Gp(x,x) which behaves as ~ (1 — x) 13 at 
x — > 1 with larger (3 than the present model Gp(x, x) = K a f a (x), as suggested by a recent 
study on p^p — ► nX [26] , one would not have got such a strongly rising behavior and would 
have obtained a negative asymmetry for (1)jv m the whole Zf region. It is quite natural to 
expect that the negative trend of the asymmetry (1)jv persists in the lower energy region, 
which is consistent with the recent HERMES data [7]. (Recall the difference between the 
present work and [7] in the sign convention of the relevant angles for the asymmetry, as 
stated just below (96).) 

It has been known that the effect of the soft gluon radiation is important in the shown 
region. However, recent study on these effects for the twist-2 double spin asymmetry based 
on the qx resummation formalism shows that the resummation effect more or less cancels 
between the polarized and unpolarized cross sections and the fixed order perturbative result 
for the asymmetry gives a reasonable prediction for the asymmetry [45]. Although the 
resummation formalism for the twist-3 cross section has not yet been developed, it should 
be the one which resums logarithmically enhanced contributions due to the emission of the 
soft gluons, whose coupling is generically blind to spin degrees of freedom. Accordingly, we 
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expect the result shown in Fig. 11 provides a correct order of magnitude for the SSA (97). 

5 Summary and conclusion 

In this paper, we have studied the SSA for the pion production in SIDIS, ep" — > enX, in the 
framework of the collinear factorization. This study is relevant for the pion production with 
large transverse momentum. In this framework, the SSA is a twist-3 observable, and the 
single-spin-dependent cross section occurs as pole contributions from internal propagators 
in the partonic subprocess, which produce the strong interaction phase necessary for SSA. 
We have derived the complete cross section formula associated with the twist-3 quark-gluon 
correlation functions for the transversely polarized nucleon, by including all types of pole 
(hard-pole, soft-fermion-pole and soft-gluon-pole) contributions. To accomplish this study, 
we have proved that the hard part from each pole contribution satisfies certain constraints 
from the Ward identities for color gauge invariance, and also clarified a complete set of the 
gauge-invariant, twist-3 correlation functions. These new developments are crucial to prove 
the factorization property for the twist-3 single-spin-dependent cross section in manifestly 
gauge-invariant form, and the "uniqueness" of the corresponding factorization formula, 
which was missing in the literature. Based on the soft-gluon-pole approximation, we have 
presented a simple estimate for the azimuthal asymmetry of SSA. It turned out that the 
non- derivative term of the SGP contribution causes sizable correction to the results using 
the derivative term only, for the small as well as moderate values of Xbj and Zf. 

Our formalism of the twist-3 calculation can be extended to study SSAs in other pro- 
cesses of current interest such as p^p — > jX, p^p — > rcX, which will be reported in the 
future publication. 
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